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Abstract. Tumour cells have to acquire a 
number of capabilities if a neoplasm is to 
become a cancer. One of these key ca- 
pabilities is increased motility which is 
needed for invasion of other tissues and 
metastasis. This paper presents a qualita- 
tive mathematical model based on game 
theory and computer simulations using 
cellular automata. With this model we 
study the circumstances under which mu- 
tations that confer increased motility to 
cells can spread through a tumour made 
of rapidly proliferating cells. The analy- 
sis suggests therapies that could help pre- 
vent the progression towards malignancy 
and invasiveness of benign tumours. 

PACS. XX.XX.XX No PACS code given 

1 Introduction 



Carcinogenesis is the process that describes the transfor- 
mation of healthy cells into cancer cells. Cancer has long 
been recognised as an evolutionary disease jT] . It is also a 
disease in which the environment (other cells in the tissue, 
nutrients, pH, etc) determines which genetic mutations 
lead to phenotypes that spread through the population. 
Hanahan and Weinberg [2] have described six capabilities 
that tumour cells have to acquire during carcinogenesis for 
a tumour to become a cancer, that is, a tumour capable of 
invasion and metastasis. These capabilities include: unlim- 
ited replicative potential, environmental independence for 
growth, evasion of apoptosis, angiogenesis and invasion. 
This paper focuses on how a tumour made of rapidly pro- 
liferating cells acquires invasiveness with some of the cells 
becoming capable of motility. We call proliferative cells 
those cells that do not rely on the environment for growth, 
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which in many cases involves mutations in genes such as 
APC, K-RAS and TGF. The motile cells are those tumour 
cells that, through abnormal regulation of the synthesis of 
cadherins and integrins, have an increase in their motility 
and invasiveness compared to other cell types [314] . This 
transition is a prerequisite for invasion and is key in the 
progression of a tumour towards malignancy. 

Tumours are the result of somatic evolution: the mu- 
tations that promote tumour initiation lead to an ecosys- 
tem in which different tumour cell phenotypes compete 
for space and resources [S]- Evolutionary game theory is 
an appropriate modelling tool in which to frame tumouri- 
genesis [6] and in which to study the possible equilibria 
between two or more phenotypes under different microen- 
vironmental circumstances. In this paper we also describe 
an agent-based model using a cellular automaton. This 
alternative model, in which the agents can be either pro- 
liferative or motile tumour cells, has been designed to mir- 
ror, as closely as possible, the game theoretical model but 
whereas the first model is suitable for analysis, the sim- 
ulations provided by the second have been used to study 
the influence of space on the results. Both models show 
that abundance of resources reduces the likelihood of suc- 
cess of motile cells. They also show that typically, there 
will be many circumstances that could lead to polyclonal 
tumours that would include both: proliferative and motile 
cells. 



2 Previous studies 

Game theory (GT), first formalised by von Neumann and 
Morgenstern [7], has a long tradition in the field of econ- 
omy [5]. It has also been successfully applied to evolu- 
tionary biology [9|10j . Its use in the field of theoretical 
medicine is much more recent although its future looks 
quite promising [BJ. GT is a powerful tool when studying 
the interactions between a number of entities called play- 
ers in which the fitness of each player depends not only 
on what the player decides to do (its strategy) but also on 
what the other players do. In evolutionary GT the strat- 
egy of an individual is an aspect of its phenotype that 
affects its behaviour in the game. One important concept 
in evolutionary game theory is that of an evolutionary 
stable strategy (ESS). A strategy is ESS if, when adopted 
by a population, it cannot be successfully invaded by an 
alternative strategy as a result of evolution. An overview 
of the application of GT in cancer research can be found 
in [11] . The first use of GT in the field of cancer research 
was described by Tomlinson and Bodmer [12)13] . Their 
models are based on pairwise interactions between dif- 
ferent types of tumour cell phenotypes in the context of 
several cancer related problems such as angiogenesis, eva- 
sion of apoptosis or the production of cytotoxic factors. 
In these problems, game theory is employed to analyse 
the circumstances that lead to polymorphism that is, co- 
existence of two phenotypes. Subsequent research carried 
out by Bach and colleagues extended this idea to interac- 
tions between three players in the angiogenesis problem 
[14] and also studied the effect of spatial dynamics on the 



equilibria [15] . Gatenby and Vincent adopted a game the- 
ory approach heavily influenced by population dynamics 
to investigate the influence of the tumour-host interface 
in colorectal carcinogenesis |16|17] and suggest therapeu- 
tic strategies [IS]. More recently Mansury and colleagues 
[19] have produced a modified game theory approach in 
which they study how the interactions between two types 
of phenotypes, proliferative and motile, in a tumour af- 
fect a number of features of the tumour growth dynamics. 
This model does not contemplate evolutionary dynamics 
and thus cannot address questions related to the effect 
of model parameters on the coexistence of different phe- 
notypes in a polymorphic tumour. Recently, Axelrod and 
colleagues have suggested, although never demonstrated, 
how GT could be used to study the possible coperation of 
different tumour sub-populations during tumour progres- 
sion [20121] . 

Cellular automata (CA) are another invention of von 
Neumann [22] who, together with Stanislaw Ulam, cre- 
ated them to study self replication. CA are based on the 
concepts of discrete space and time and are well suited to 
model processes at the cellular level [23]. CA have been 
extensively used to model many aspects relevant in can- 
cer research |24|25] such as early tumour growth [19|26j . 
invasion[27 28 29 30 , angiogenesis [31132] and evolution 
[33134] , 

3 The game theory model 

The model assumes a very simple microenvironment in- 
habited initially only by proliferative cells mimicking the 
situation in many benign tumours [1]. The motile pheno- 
type can appear as a result of one or more genetic mu- 
tations. The interactions between these two phenotypes 
are characterised in table [T] Parameter b is the base pay- 
off that corresponds to a cell that does not need to share 
resources (glucose, oxygen, etc) with any one else. The 
cost of motility c represents the fitness loss incurred by 
any motile cell that at a given moment chooses to move 
to a neighbouring location. This cost could represent for 
instance the fact that cells that are capable of motility 
cannot divide while moving [35136] . When a motile cell 
finds a proliferative cell then the proliferative cell stays 
and uses all the resources of that location. The motile cell 
moves to a different location where it is assumed that it 
will find a similar amount of resources. When two motile 
cells meet then one of them stays and takes the resources 
whereas the other moves to another location. When two 
proliferative cells interact they have to share the resources 
of the location in which they find themselves. Given the 
payoffs shown in tableHJ the Bishop-Cannings theorem [9] 
can be used to obtain the mixed strategy that would be 
an evolutionary stable strategy (ESS). A mixed strategy 
is a strategy in which one of a number of pure strategies 
can be used with a given probability. If the pure strate- 
gies are to be motile or to be proliferative then a mixed 
strategy could be 0.6 implying that the cell would play the 
strategy of motility with probability 0.6 and the strategy 
proliferative with 0.4. A mixed strategy in this game could 
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Fig. 1. The graph shows the results of analysing the game 
theoretical model. The X axis shows the value of the cost of 
motility and the Y axis shows the proportion of the population 
with the motile phenotype. 

also correspond to a population composed of more than 
one phenotype. 



Table 1. Payoff table that represents the change in fitness of 
a tumour cell with a given phenotype interacting with another 
cell. The base payoff in a given interaction is b and the cost of 
moving to another location with respect to the base payoff is 
c. 
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The genetic polymorphism, in which the two pheno- 
types coexist in the tumour, can be calculated if we as- 
sume that for an equilibrium to exist, the payoff of one 
strategy is equal to the payoff of the other strategy, thus 
pE(m,m) + (1 —p)E(m,p) = pE(p,m) + (1 — p)E(p,p), 
where p is the probability of a cell having the motile phe- 
notype and E(x,y) is the payoff for a cell with phenotype 
x when meeting a cell with phenotype y. The values of E 
can be extracted from table [1] and thus, the probability 
p of adopting strategy motility in equilibrium varies with 
b and c (with c not more than half the base payoff) as 
shown in the following equation: 



Fig. [T] shows the proportion of motile phenotypes in 
the population when the base payoff, b, is fixed. It can be 
proven that in games with only two pure strategies there 
is always one stable mixed strategy and the composition 
of pure strategies in it corresponds to a polymorphism [9] . 



4 The cellular automaton model 

The GT model does not take into account the influence of 
space in the dynamics of tumour evolution. For that rea- 
son we have introduced a cellular automaton (CA) with 
the aim of producing a spatial model in which the be- 
haviour of proliferative and motile cells can be investi- 
gated. Space in the CA is characterised by a 1000 x 1000 
lattice in which each site can host one or more tumour 
cells. The carrying capacity of a lattice site corresponds 
to the area in which one cell can affect another one. Each 
lattice site also contains nutrients which are needed for cell 
proliferation. Nutrients are assumed to be homogeneously 
distributed in the lattice and replenished to a given con- 
stant value each time step. Thus the resources that deter- 
mine cell proliferation are available space and nutrients. 
Cells interact and share resources only with other cells 
in the same lattice site. This very simple microenviron- 
ment mimics the one that is implictly assumed in the GT 
model. The behaviour of the tumour cell types was de- 
signed to resemble, as closely as possible, that of the two 
phenotypes of the GT model. It features proliferative and 
motile cells. The former do not move even if that means 
sharing an ever decreasing amount of resources. The latter 
divide only when they do not have to share the available 
resources with no other cell and try to move to a better 
location otherwise. The precise behaviour of the different 
phenotypes is the following. The cells with the prolifer- 
ative phenotype divide if there are enough nutrients in 
the lattice site. The new offspring will be created in the 
same site as the progenitor unless the population of cells 
at that site has exceeded the maximum carrying capacity. 
In the last case the offspring is created in a neighbouring 
location having the lowest density of cells (see figure [2|) . 
Motile cells will never share a lattice site with another cell 
if there is more available space in a neighbouring location. 
Thus, motile cells always move to a neighbouring site if 
the density of the new location is lower than that of the 
current location. If more than one neighbouring location 
is equally suitable then the new location will be selected 
randomly among them. On the other hand if the location 
is suitable and there are enough nutrients then the cell 
will produce offspring in the same location as the parent 
(see figure [3]). 

The update of the CA is as follows: 

— If cell is motile and there are neighbours: 

— Look for the neighbouring lattice site with the small- 
est cell count. 

— Move to new location if cell count is smaller than 
in current location. 

— Otherwise: 

— If there is enough space in the current lattice site 
then divide. 

— else 

• Look for the lattice site with the smallest num- 
ber of cells. 

• If chosen site has enough space produce off- 
spring. 

• If offspring has been produced then, with a 
small probability, mutate new cell to motile. 
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Fig. 2. Behaviour of proliferative cells, (a) In the example each 
lattice site can host up to 8 tumour cells. When the lattice site 
contains enough space for the proliferative cells to produce 
their offspring in their current location they will do so. (b) If 
the lattice site can not contain the new offspring then it will be 
created in the neighbouring site with the most available space. 
In the example the first two cells selected have enough space in 
their site to create their offspring. The remaining ones have to 
select a neighbouring site and since empty ones are preferred 
then each cell choses a different neighbouring site. 



(a) ts=0 



(b) ts=100 



(c) ts=200 





(d) ts=300 



(e) ts=400 



(f) ts=500 



Fig. 4. Result of a sample simulation produced by the CA 
model. In this particular simulation each lattice site can host 
up to 10 tumour cells. Every time step only one unit of nu- 
trients is replenished in each site. These circumstances reduce 
the cost of motility and thus select for the motile phenotype. 
Lattice sites which contain exclusively proliferative cells are 
marked in blue, those that contain only motile cells are red 
and finally those that contain more than one type of cell are 
orange. 



(a) (b) 

Fig. 3. Behaviour of motile cells, (a) In this example motile 
cells are red and proliferative cells are blue. If a motile cell is 
not alone in the lattice site then it will move to a neighbouring 
site with the smallest cell density. This will leave all the site 
resources to the proliferative cell, (b) On the other hand if it 
does not have to share the site with any other tumour cell then 
it will divide. In the following time step, at least one of them 
will have to move to a different location if there is any site with 
a lower cell count in the neighbourhood. 

Every time step each of the tumour cells is chosen to 
have its programme executed. The order in which the cells 
are called is randomly selected in each time step to avoid 
potential biases. The simulations start with a group of 
proliferative cells placed in the centre of the lattice (cir- 
cumference with radius 20 lattice sites). Every time pro- 
liferative cells divide there is a small chance (10~ 5 ) that 
there will be a mutation conferring the offspring the motile 
phenotype. Motile cells on the other hand will never pro- 
duce proliferative cells when dividing. The reasoning for 
this is that the mutations in our model are not necessarily 
genetical but all have an effect on the phenotype. Whereas 
cells in tumour in an early stage of development tend to be 
proliferative, the acquisition of motility and invasiveness 
is a necessary capability that may emerge as a result of a 
number of genetic and epigenetic mutations [2J. Hence, it 
is very unlikely that a motile cell could produce offspring 



that would accumulate the necessary mutations to disable 
this capability. 



5 Simulations and analysis 

The CA model was used to produce several simulations 
(e.g.: figured]), each one lasting 5000 timesteps, in which 
different values of nutrient replenishment and capacity of 
the lattice sites were tested. 

Figure [5] shows the results of the simulations. The 
curves represent the proportion of motile cells that re- 
sults for values of nutrient replenishment between I and 
25. Each line represents a different lattice site carrying 
capacity of 10, 20, 50 and 100 cells per site. 

The results from both analysis and simulations support 
current medical knowledge: if the cost of motility is zero 
then the motility phenotype will reach fixation and spread 
through the tumour population. If the cost of motility in- 
creases, so does the proportion of proliferative cells in the 
population. The proportion of motile cells in the simula- 
tions depends on the carrying capacity and the amount of 
nutrients. Simulations in which lattice sites can host many 
cells but nutrient availability is scarce, the cost of motility 
in relation to the cost of remaining in the same place is 
quite low and thus yield a high proportion of motile cells. 
On the other hand when lattice sites can carry just a few 
cells forcing proliferative cells to produce their offspring 
in neighbouring sites then the relative cost of motility in- 
creases and likewise the proportion of proliferative cells. 
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Fig. 5. The graph shows the proportion of motile cells (Y 
axis) that result from using different amount of nutrients per 
time step (X axis) after running the CA model for 500 time 
steps. The results are for different carrying capacities of each 
lattice site (10, 20, 50 and 100 cells/site). In all cases increasing 
the amount of nutrients supplied each time step reduces the 
proportion of motile tumour cells. 



In the CA model the ratio of nutrients per lattice site 
capacity gives a measure of the cost of motility (c) rela- 
tive to the maximum fitness benefit (b) . When the ratio of 
nutrients per lattice site capacity is high then the relative 
cost of motility will be high too. In these situations there 
are plenty of nutrients so that proliferative cells normally 
have enough resources to divide in every time step and 
thus out-compete motile cells. On the other hand if the 

ratio of nutrients per cell site is low then many prolifera- Acknowledgements 
tive cells will not have enough nutrients and might have 
less opportunities to produce offspring than cells with the 
motile phenotype. 



and nutrients) should have an influence of the fitness of 
the two phenotypes. Many of these simplifications are nec- 
essary and similar to those adopted by other researchers 
!12|13|14|15TT9] . This fact does not detract merit from the 
potential of the model to test the ideas behind the word 
models usually put forward by oncologists [2]. The simu- 
lations also support a therapy that would raise the cost 
of motility with respect to the base benefit to reduce the 
proportion of tumour cells capable of motility. Thus, by 
increasing the relative fitness advantage of non-motile phe- 
notypes increases the probability of the tumour remaining 
benign. Increasing the cost of motility could be achieved 
in at least two different ways. One option would be to de- 
sign therapies to hinder the detachment of tumour cells 
from the extra-cellular matrix by means of downregula- 
tion of integrins [37]. Alternatively, the model suggests 
the counter-intuitive idea that it should be possible to 
increase the fitness of proliferative cells over motile cells 
by increasing the amount of nutrients in the neighbour- 
hood of the tumour. Increasing the amount of nutrients 
(equ. QJ produces the same effect as acting directly on the 
cost of motility. It seems reasonable to think that increas- 
ing the amount of nutrients would act as a disincentive 
for a tumour cell to become motile. This is not neces- 
sarily a therapy that would work in all cases: there are 
benign tumours that are life threatening even if they do 
not become invasive. Despite that, there are many cases 
in which a growing but non aggressive tumour will have a 
much better prognosis than a smaller but invasive one. In 
those cases this therapy could help the patient by produc- 
ing a tumour environment that selects for less aggressive 
phenotypes. 



6 Conclusions 

We have introduced a game theoretical model that allows 
the study of the emergence of invasion and motility in 
tumours made of proliferative cells. We have also shown 
a CA model designed to study how spatial considerations 
could alter the results yielded by the GT analysis. The ex- 
istence of polymorphic tumours has been found in other 
GT models [13] and the results from both, our GT analy- 
sis and the CA-based simulations, confirm that there are 
many circumstances under which different strategies could 
be present in the tumour cell population at a given time. 
Many simplifications were made in order to get a model 
tractable enough to be simulated and analysed. It is clear 
that more than two phenotypes are possible at this stage 
of carcinogenesis and that tumour cells usually interact in 
a more complex fashion than the one used in this study. 
The microenvironment was also kept as simple as possi- 
ble. It is clear that more complex microenvironments (with 
differences in the concentration of e.g. pH, oxygen, ECM 
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